Orbital-resolved Soft X-Ray Spectroscopy in NaV 2 C>5 
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We demonstrate that angle-resolved soft x-ray spectroscopy can resolve absorption by inequivalent 
oxygen sites and by different orbitals belonging to the same site in NaV^Os. By rotating the 
polarization direction, we see a dramatic change in the absorption spectra at the oxygen K edge. 
Our theory identifies the detailed composition of the spectra and predicts a correct energy-ordering 
of the orbitals of three inequivalent oxygen atoms. Because different orbitals dominate absorption 
spectra at different energies and angles, one can excite at a specific site and "orbital". In contrast, 
absorption at the vanadium L edge does not show large changes when varying the polarization 
direction. The reason for this is that different excitation channels (involving different initial states 
for the excited electron) overlap in energy and vary in compensating ways, obscuring each channel's 
sensitive polarization dependence. 

PACS: 78.70.Dm, 71.15.-m 



I. INTRODUCTION 



Orbital degrees of freedom play a crucial role in many technologically interesting materials, such as superconductors, 
magnetic materials, and other strongly correlated compounds. For instance, the spin-orbit coupling effects determine 
whether the magnetic crystalline anisotropy in magnetic thin films favors magnetic moments normal or parallel to the 
surface of the sample. Suchdmportant effects motivate many investigations that probe orbitals in real materials. Soft 
x-ray spectroscopy (SXS)uij has a unique capability to reveal spin and orbital information specific to one element. 
By employing x-ray magnetic circular dichroism, one can measure spin and orbital moraentsp using x-ray magnetic 
linear dichroism, one can image antiferromagnetic domains at sub-micrometer scales.13 Despite these tremendous 
achievements, imaging atomic orbitals has not often been done. The first attempt at the C K edge of graphite 
done by Rosenberg et alu A more elaborate study was performed a decade later with an intense synchrotron source 
On surfaces, the first investigation was made by Nilsson et aZ.II2l They used s-polarized soft x-rays to selectively excite 
7r orbitals of CO adsorbed on Ni and used p-polarized light at off-normal angles to excite both a and 7r orbitals, 
from which the orbital information could be extracted. However, obtaining such valuable orbital information in bulk 
materials is still a challenge for SXS. 

In this work, we consider NaV20s, a fascinating transition- metal compound that features three inequivalent oxygen 
sites: Oi, O2 and O3 [cf Fig, r-Ua)]. At 34 K, this material undergoes a spin-Peierls-like phase transition, which 
has attracted great interest .Ej~L3 NaV^Os is orthorhombic with symmetry of P mm n and lattice constants a=11.316 
A, 6=3.611 A and c=4.797 A. It is a ladder structure compound along the b axis with V-Oi-V bonds as its rungs 
extending along the a axis. The neighboring ladders along the b axis arc connected by shifting a half-unit cell along 
the b axis. These ladders form quasi-two dimensional layers between which the sodium atoms are intercalated. One 
Oi, three oxygen O2 and one apical O3 atoms form a pyramid inside which one vanadium atom is situated. These 
three different oxygen atoms present a challenge to site-specific SXS, because one desires to obtain information, not 
only about different orbitals of atoms of different elements, but also about different orbitals on the same site. 
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Specifically, we report orbital-resolved x-ray absorption measurements, using the total electron yield method (TEY), 
on NaV2C>5 at the O K (Is) and V L3 (2p 3 / 2 ) edges. By changing the angle 9 between the electric field-polarization 
direction and certain crystal axes of NaA^Os, we can determine characteristic orbital energies, not only for orbitals 
on inequivalent O sites, but also for different orbitals on the same O site. With the aid of theoretical calculations, 
we are able to pinpoint the origin of the angular dependence of the O K edge absorption spectrum, distinguish 
different contributions thereto, and accurately determine the crucial energy ordering of the different orbitals from 
three inequivalent oxygen atoms. Combining different scans along different axes, we are able to project out spatial 
orientations of oxygen orbitals-in real space. Thus, our study highlights an important capability of SXS that has not 
been explored in great detail.li3 

Our experiment reveals a surprising difference between the vanadium and oxygen edges. At the V L3 edge, the 
spectra show a very small dependence on 9, but at the O K edge there is a large dependence. Our calculations are 
scalar-relativistic, so that we model the V £2,3 manifold without spin-orbit effects (including P1/2—P3/2 splitting) taken 
into account. However, this description should still be able to address many aspects of a polarization dependence, 
which depends on states near the Fermi level. The theory shows that, subject to the dipole selection rule Al = ±1 
and Am = 0, ±1, transitions from s orbitals only have one initial state, i.e., m — 0, while transitions from p orbitals 
have three initial states, i.e., m = 0, ±1. Polarization dependences for absorption by electrons in each of the initial 
states compensate each other and ultimately lead to a very small polarization dependence. This explains the observed 
trend for the V L3 (2p 3 / 2 ) edge. 

The remainder of this article is arranged as follows. A brief review of our experimental measurements is described 
in Section II. In Section III, we present our theoretical scheme. The results and discussions are shown in Sec. IV. 
Finally, we conclude our paper in Section V. In a companion paper ,E£I we analyze emission spectra and experimentally 
determine the occupied p density of states from each oxygen site. This paper is devoted entirely to the analysis of 
absorption spectra. 

II. EXPERIMENTAL MEASUREMENTS 

Single crystals of NaV2 05 were prepared at the Oak Ridge National Laboratory.!!!! The crystals grow as rectangular 
platelets with the c axis normal to the plate surface, the a axis along the short side of the rectangle and b along the 
long side, as shown in Fig. 1(b). The samples were cleaved along the c axis. The experimental measurements were 
done at Beamline 8.0 of the Advanced Light Source at the Ernest Orlando Lawrence Berkeley National Laboratory.lla 
The geometry of the measurements is indicated in Fig. 1(b). Horizontally polarized x-rays from an undulator and 
monochromator were incident on the sample. In the TEY detection mode, the x-ray absorption spectrum is measured 
by collecting the electrons that escape from the sample. Two samples were used in the experiment, and mounted at 
right angles to each other. Rotation of the sample holders scanned the polarization vector of the incident x-rays in 
the a-c plane for one sample and in b-c plane for the second sample. 

NaV205 is an ideal system, in which to reveal a polarization dependence, because its structure is highly anisotropic. 
In Fig. 2, we show the normalized absorption spectra of a polarization scan from along the b axis to along the c axis, 
specifically from 9 = 7.5° to 9 — 75° with a step of 7.5°. For each angle, the normalization is done with respect to 
the total area between 514 eV and 536 eV because this energy region covers both the V-L and O-K edges. 

We note that the TEY method has a finite sampling depth of about 100 A. This suggests that our electron yield 
arises from a local spot of the sample. To illuminate the sample in the same place during measurement and to ensure 
data consistency, we adjusted the sample position after each change of its orientation angle. We found that data 
obtained using this technique are reproducible within 5 %. 

We first consider the vanadium L3 edge around 518 eV [Fig. 2(a)]. When we scan the polarization from along the 
b axis to along the c axis, the spectra change very little. We observed a similar behavior when we scanned from along 
the a axis to along the c axis (not shown). However, for the oxygen K edge, we see a totally different picture [Fig. 
2 (b)]. Increasing the angle 9 from 7.5° to 75° strongly changes the spectra. Three main peaks A, B and C, located 
at 530.6 eV, 532.1 eV and 532.8 eV, respectively, can be clearly identified. Peak A's intensity decreases sharply and 
eventually disappears as the polarization moves to along the c axis, which is normal to the sample surface. In the 
meantime, intensities of peaks B and C increase strongly and become the two main peaks on the spectra. To have 
a clear view of the intensity change, in Fig. 3 we plot the intensities versus angle 9 for those three peaks. Their 
variations can be nicely fitted to a cos 29 dependence, which is the simplest possible variation (see the solid lines in 
Fig. 3). 
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III. THEORETICAL FRAMEWORK 



Our theoretical calculation to simulate the experiment is based on pseudopotentialJ3 local-density-approximation 
(LDA) calculations Jla Electron states are found by solving the self-consistent Kohn-Sham (KS) equations, 

h 2 

ff*nk(r) = {- — V 2 + Kff }*„k(r) = £U*nk(r) (1) 
2m 

Here \l/„k is the Bloch state for an electron in band n at crystal wave-vector k. The effective potential is given by 

VeS = \^on + Vn(r)+V xc (r). (2) 

The three terms account for the pseudopotential interaction of electrons with ion cores Vi on , the Hartree potential 
Vh(i") that is obtained from the electron charge density, given by 

occ'd. 

n(r) = l*«k(r)| 2 (3) 

nk 

and exchange and correlation effects V xc . We use the Ceperley- Alder formEl for the exchange term V xc — 
SE xc [n]/5n(r). The ionic potential Vion is a norm-conserving pseudopotential computed within the HartreepEock 
approximation, and with core-valence correlation effects included using the core-polarization-potential method.E3 As 
basis sets, we use plane waves or optimized basis functions^ _ — . _ — . 
To compute an x-ray absorption spectrum, we use the expression, EiO 

SH(x27r^|</|p-A|i>| 2 S(u + Ei-E f ) (4) 
/ 

where p • A is the electron-photon interaction approximated by a dipole operator; \i) is the initial core state with 
energy Ei, /} is the final state with energy Ef, and u> is the incident photon energy. In our present analysis, which 
mainly deals with determining the symmetries of orbitals, we neglect electron-core hole interactions in the final state. 



IV. RESULTS AND DISCUSSIONS 



A. Oxygen K (Is) edge absorption spectrum 

In the experiment, we always rotated our polarization direction from along the b axis to along the c axis in the 
way already described. Theoretically, we did calculations in the same way. In Fig. 4(a), we show theoretical x-ray 
absorption spectra at the O K edge, where three peaks, A, B and C, are also clearly observed with different intensities 
as a function of polarization. The spectra were calculated with 12 x 12 x 12 k-point meshes. The results are found 
well-converged. The angular dependence of these peaks is in good agreement with the experimental observations: 
Peak A decreases as the polarization moves away from along the b axis, whereas both peaks B and C increase. The 
change also follows a cosine function of 9, consistent with the experimental observation. Note that the experimental 
measurements cannot directly determine which oxygen orbitals contribute to these main peaks, whereas calculations 
should provide important insight into this. 

Figures 4(b-d) show the respective contributions from three inequivalent oxygen sites. To simplify our notation, we 
let O m (P ra ) denote the nth orbital of oxygen O m that participates the absorption process. That is, the local, partial 
density of states for such an orbital can play a significant role in absorption as a function of excitation energy, and 
this role can be strongest at some "characteristic" energy. We first consider peaks B and C. By comparing Fig. 4(a) 
with Figs. 4(b), 4(c) and 4(d), we infer that peaks B and C in Fig. 4(a) and Fig. 2(b) originate from Oi(Ps), 02(P4), 
(^(Ps), and 03^5) orbitals, whose characteristic orbital energies (expressed relative to the O Is level) are 533.0 eV, 
533.1 eV, 533.7 eV and 533.1 eV, respectively. 

When considering peak A in Figs. 4(a) and 2(b), we notice that all three inequivalent atoms contribute but are 
favored in different ranges of energy. Figure 4(b) illustrates the main contributions from Oi. It is important to note 
that only Oi(Pi) appears at 530.2 eV, which helps us to resolve different orbital contributions (see next paragraph). 
The peak intensity decreases as the polarization direction deviates from the b axis, in a similar way as peak A in Fig. 
2(b) or Fig. 4(a). Oi(Pi) is immediately followed by C>3(Pi) at 530.6 eV with an energy difference between them 
of about 0.4 eV, [cf. Fig. 4 (d)]. The Oa(Pi) intensity changes in a similar manner as does that for Oi(Pi). After 
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the C>3(Pi) orbital, the 02^3) orbital gradually appears around 531 eV to 532 eV [cf. Fig. 4(c)]. Because of their 
unique positions, O2 orbitals are distorted, and not along one particular axis. Thus, the angular dependence of their 
contribution to absorption is weaker than for orbitals on Oi and O3. 

These results show that the measured peak A has contributions from the Oi(Pi), 03(Pi), and part of 02^3) 
orbitals. This could prove crucial for our tentative explanation of emission spectra that is given subsequently in Ref. 
16. In particular, because our theory shows that these three orbitals have different energies, it is possible to excite 
different orbitals at different energies and to resolve them experimentally. This is exactly what we did. Because the 
energy increases in the order of Oi, O3 and O2, we can selectively excite these atoms. When we excite at the Oi(Pi) 
orbital edge, we only see one peak in the emission spectrum, because O2 and O3 atoms are not excited, though in the 
normal fluorescence spectrum, we expect to see two main peaks. Therefore, with theoretical guidance, we demonstrate 
that we are able to resolve these different orbitals associated not only with different sites such as Oi(Pi), 3 (Pi), 
02^3) but also with the same site such as Oi(Pi) and Oi(P3). 

By scanning polarization from along the a axis to along the b axis and from along the a axis to along the c axis, we 
can project out orbital orientations in real space. Putting together all the results, we can show some of the orbitals in 
the real space in Fig. 5. One sees that there are two kinds of O3 orbitals: one is along the a axis, the other is along the 
c axis, and one Oi orbital along the b axis, two O2 orbitals along the b and a axes. Note that, unlike charge-density 
plots for orbitals, the orbitals presented here are experimentally accessible. Thus, we demonstrate that the soft x-ray 
spectroscopy has the capability of resolving orbitals from the same site. 

B. Vanadium L s (2p ;j / 2 ) edge 

The x-ray absorption spectrum at the V L3 edge is rather unique and shows a very weak angular dependence. At 
room temperature, there is only one inequivalcnt vanadium atom. In our scalar-relativistic calculations, excitation 
of V 2p electrons can involve three electron initial states: m = 0, ±1. Our absorption spectrum intrinsically involves 
a sum over excitations with the indicated initial states, and the experimentally observed spectrum is an analogous 
sum at the V L 3 edge. We first show the results of the first channel in Fig. 6(a). We note that with the incident 
polarization turning away from the b axis, the whole spectrum clearly changes. If this were the dominant excitation 
channel, then one would see a strong angular dependence in experiment. However, our experiment showed otherwise, 
which implies that we should take into account the other channels. In Fig. 6(b), we plot the change because of the 
I = 1, m = ±1 channel. Without spin-orbit coupling, the spectra for both I = l,m = +1 and I = 1, m = —1 channels 
are same. Two prominent features are: (1) the mean peak strongly overlaps the spectra of the I = 1, m = channel; 
(2) the portion showing strong angular dependence is on the right-hand side, which is opposite to that in Fig. 6(a). 
Because of this compensating effect, the total absorption spectrum does not show a strong angular dependence as one 
can see in Fig. 6(c). This explains our experimental observations. 

Now comparing this with the oxygen K-edge spectra where only one excitation channel, I = and m = 0, 
participates in the process, we clearly see that appropriate summation over the multiple excitation channels (i.e., 
corresponding to different electron initial states) obscures the polarization dependence, which ultimately leads to the 
big difference between the role of polarization at the V L3 edge vs the O K edge. This finding demonstrates the 
importance of understanding selection rules in x-ray spectroscopy, and it might help one optimize the polarization 
direction to reveal the maximal information from x-ray spectroscopy measurements in the future. In particular, it 
would be interesting to ask whether the present polarization effects for the various components of the V L3 spectrum 
could be resolved by circular dichroism. However, the difficulty is that even with the 2p core spin-orbit coupling, in 
the absence of spin-polarization of the V-d bands, Am = ±1 transitions should also be the same. 

V. CONCLUSIONS 

In conclusion, we measured angular dependence of x-ray absorption spectra at the vanadium L3 edge and oxygen 
K edge. Rotating the polarization from along the b axis to along the c axis, the spectra at these two edges behave 
quite differently There is very little change at the V L3 edge, while a dramatic change is observed at the O K edge. 
Such a large change enables us to selectively excite different orbitals belonging to the same site. Thus, we are able to 
resolve orbitals for the same element using soft x-ray spectroscopy. Our theoretical calculation provides the detailed 
information to understand the experimental spectra. It clarifies how the predominant first peak at the oxygen K 
edge is from three different orbitals of three different oxygen sites (Oi, O2 and O3) and determines the crucial energy 
ordering of the orbitals on these sites. This also helps explain the emission spectra. Our calculations also show how 
a compensating effect can arise in the polarization dependence of various excitation channels for the V L 3 edge, so 
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that the total absorption spectrum at [the edges exhibits a small polarization dependence. After we submitted this 
paper, we noticed three recent papersEj that intensively debated whether the anomalous scattering factors in the 
x-ray diffraction technique could be used to probe the charge ordering in NaV2 05. Our study shows that if excited 
at the V L3 edge, the orbital or structural effects are small. Thus, one might explain the anomalous scattering as 
resulting from the charge ordering, provided that the scattering intensity of the charge ordering is larger than that of 
the orbital ordering. 
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FIG. 1. (a) NaV20s structure; (b) Experimental geometry. 



FIG. 2. Experimental x-ray absorption spectra of NaV20s at (a) the V L3 edge and (b) the O K edge. The angle 6 is 
between the electric polarization and the b axis. 6 increases from 7.5° to 75° in steps of 7.5°. The scan is in the b-c plane. The 
vertical dashed lines denote the positions of peaks A, B and C at 530.6, 532.1 and 532.8 eV, respectively. 



FIG. 3. Intensities versus the polarization angle 9 for peaks A, B and C in Fig. 2(b). 

FIG. 4. Theoretical absorption spectra of NaV20s at the O K edge, (a) Total absorption spectra at the O K edge. 
Contributions are from (b) Oi, (c) O2 and (d) O3. The scan condition is same as in Fig. 2. The arrows denote the change in 
intensity while varying polarization from along the b axis to along the c axis. 



FIG. 5. Oxygen orbitals' spatial distributions in NaV20s. 

FIG. 6. Theoretical V L3 edge absorption, (a) Only the channels / = 1 and m = is calculated; (b) Only the channel / = 1 
and m = ±1 are calculated, (c) Total absorption spectrum at the V L 3 edge. The scan condition is exactly same as in Fig. 2. 
The arrows have the same meaning as in Fig. 4. 
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